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ABSTRACT 
A  forecast  model  proposed  by  'Arnason  as  being  capable  of  producing 
realistic  forecasts  of  displacement  of  systems  at  and  below  500  mD  j_s 
investigated.  The  prediction  equation, 


o 


is  applied  to  test  cases  at  1000,  850,  700  and  500  mD  with  various  values 
of  the  parameter  k]_.  The  effect  of  the  value  of  the  parameter  on  the 
amount  of  divergence  present  in  the  model  is  discussed.  Root-mean-square 
errors  for  2li-hour  and  l|8-hour  forecasts  with  various  values  of  k-^  at 
each  level  are  presented.   Optimum  values  of  k->  for  the  cases  tested  are 
selected  as  3.5,  2.75  and  0.9  for  the  levels  850,  700  and  500  mb  respec- 
tively. It  is  found  that  with  proper  values  for  k-,  the  model  will  provide 
good  forecasts  of  system  movement  at  each  level  tested,  but  as  expected, 
changes  in  the  intensity  of  systems  are  not  correctly  forecast.  In  all 
cases  values  of  root-mean- square  error  of  forecasts  are  less  than  those 
obtained  by  applying  the  "Helmholtz"  barotropic  model.  A  geostrophic 
version  of  the  model  is  also  tested  and  is  found  to  produce  results  simi- 
lar to  those  obtained  for  the  stream-function  version  except  that  the 
geostrophic  "blow-up"  of  low- latitude  high  pressure  cells  is  observed. 

The  writer  wishes  to  express  his  appreciation  for  the  guidance  and 
encouragement  given  him  by  Professor  George  J.  Haltiner  of  the  U.  S. 
Naval  Postgraduate  School  in  the  course  of  this  investigation.  In  addi- 
tion, appreciation  is  expressed  to  Mr.  Geirmundur  'Arnason,  on  whose  work 
this  paper  is  based,  and  to  personnel  of  the  Fleet  Numerical  Weather 
Facility  for  many  instances  of  assistance. 
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1.   Introduction 

Since  the  initation  of  operational  numerical  weather  forecasts  by 
the  Joint  Numerical  Weather  Prediction  Unit  in  1955?  the  only  forecast 
model  which  has  stood  the  test  of  daily  operational  use  has  been  the 
barotropic  model  first  proposed  by  Rossby  in  1939»  The  barotropic  model 
remained  in  its  original  form  until  1958  when  Wolff  [l]  discovered  large 
systematic  errors  due  to  improper  forecast  of  very  long  waves.  After 
these  errors  had  been  empirically  corrected  by  Wolff,  it  was  discovered 
by  Cressman  [2J  that  these  long-wave  errors  could  be  reduced  by  the  in- 
troduction of  divergence  into  the  barotropic  model.  This  resulted  in 
the  divergent  or  "Helmholtz"  version  of  the  barotropic  model  which  is  in 
general  use  at  this  time. 

With  the  improved  barotropic  model,  high  quality  forecasts  have  been 
available  for  the  500-mb  level.  This  barotropic  model,  however,  is  by 
nature  restricted  in  its  use  to  the  500-mb  level  and,  in  addition,  is  ca- 
pable mainly  of  forecasting  the  movement  of  systems,  and  not  in  general 
changes  in  intensity.  Occasionally  attempts  have  been  made  to  replace 
the  barotropic  model  in  operational  use  with  some  type  of  baroclinic 
model  which  would  obviate  these  shortcomings.  These  attempts  have  so  far 
met  with  failure  largely  due  to  the  inability  to  control  the  forecast  of 
cyclogenesis. 

•Arnason  has  suggested  that,  in  view  of  the  great  utility  of  the 
500-mb  forecasts  produced  by  the  barotropic  model,  it  might  be  worth- 
while temporarily  to  set  aside  the  problem  of  predicting  cyclogenesis 
and  to  try  to  develop  a  model  which  would  possess  the  same  skill  of  fore- 
casting movement  at  other  levels  which  the  barotropic  has  at  500  mb.  To 
this  end,  'Arnason  has  proposed  a  model  £3]  which  maintains  many  of  the 
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characteristics  of  the  barotropic,  but  which  is  more  adaptable  to  other 
levels. 

This  model  is  based  on  a  single-layer  compressible  fluid  with  a 
stratified  density  distribution,  bounded  by  two  constant -pressure  sur- 
faces. The  theoretical  development  results  in  an  expression  for  the 
horizontal  wind  divergence 

Jw\v- x.( k£  +  v«v*\  (i) 


Introducing  this  expression  into  the  vorticity  equation, 

jj  +  VV7-7  dw  V  =  0  j  (2) 

J*        ' 

gives  the  prediction  equation  for  the  "stratified"  model 


TT  /  F7TU*  y 


(3) 


In  order  to  apply  this  equation  to  forecasting,  'Arnason  has  suggested, 
as  one  alternative,  the  use  of  a  stream  function  for  computation  of  wind 
and  vorticity  and  the  elimination  of  $>  through  a  geostrophic  assump- 
tion. This  gives 


I  v*-£2  M£  nV'Vy  -£2.  V  *  V  V  =  0, 


which  is  of  the  same  form  as  the  prediction  equation  for  the  divergent 
or  "Helmholtz"  barotropic  model,  but  with  an  additional  divergence  term 
involving  V  and  with  the  substitution  of  the  stability-dependent  term 
g1  for  g. 

Although  the  forms  of  the  prediction  equations  are  similar,  there 
are  several  basic  differences  between  the  "stratified"  and  other  barotrop* 
ic  models.  By  the  method  of  perturbation  analysis  'Arnason  has  obtained 
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the  meteorologically  significant  root  of  the  frequency  equation  for  the 
stratified  model 

Although  this  equation  is  again  similar  to  the  corresponding  equation 
for  the  Helmholtz  barotropic  model,  an  important  difference  arises  in 
the  fact  that  the  basic  current  L7   does  not  have  a  multiplying  factor. 
By  empirically  varying  the  value  of  g'H,  the  magnitude  of  the  6  term 
may  be  changed  without  reducing  the  effect  of  the  basic  current,  thus 
allowing  control  over  the  speed  of  movement  of  systems.  As  can  be  seen 
from  equation  (l)  this  is  equivalent  to  a  variation  in  the  amount  of 
divergence  present  in  the  model.  This  provides  a  method  of  adjusting 
the  model  to  yield  proper  system  movement  at  a  number  of  levels. 

The  stratified  model  as  proposed  by  'Arnason  seems  to  provide  a 
promising  approach  to  numerical  forecasting  at  and  below  the  5>00-mb 
level.  It  is  the  purpose  of  this  paper  to  apply  the  model  to  a  number 
of  test  cases  and  determine  its  effectiveness. 


2.   Application  of  the  Model  to  Forecasting 

In  the  application  of  the  prediction  equation  (3)  to  actual  fore- 
casts there  are  several  choices  to  be  made.  The  parameter  V   has  not 
been  explicitly  defined,  but  may  be  taken  to  be  any  kind  of  space- 
averaged  or  zonal  wind.  Also  the  value  of  g'H  is  left  to  be  determined 
empirically. 

By  examination  of  the  prediction  equation  (3)>  it  can  be  seen 
that  the  additional  "divergence"  term  of  the  stratified  model  which  in- 
volves V  essentially  represents  the  horizontal  advection  of  height 
by  the  mean  wind.  This  immediately  suggests  the  notion  of  a  steering 
current,  which  may  be  thought  of  as  a  wind  at  some  level  or  as  a  mean 
wind  in  some  layer  of  the  atmosphere.  Also,  the  magnitude  of  the  di- 
vergence term  is  dependent  on  the  difference  between  the  mean  wind  and 
the  wind  at  the  forecast  level.  Thus  it  appears  reasonable  to  use  as 
a  mean  wind  something  close  to  the  500-mb  wind,  since  this  would  mini- 
mize the  divergence  term  at  500  mb  and  allow  it  to  increase  with  dis- 
tance from  500  mb.  This  is  generally  consistent  both  with  observed 
patterns  of  divergence  in  the  atmosphere  and  with  the  need  for  increas- 
ing the  value  of  the  divergence  term  at  lower  levels  to  supplement  the 
movement  due  to  vorticity  advection  alone. 

A  second  parameter  which  must  be  empirically  interpreted  is  the 
factor  g'H.  The  value  of  g'  is  a  measure  of  the  static  stability  of 
the  basic  flow  and  H  is  the  thickness  of  the  layer  considered  in  the 
model.  Since  neither  is  uniquely  defined,  the  combination  may  be  taken 
as  any  reasonable  value  which  produces  proper  movement  of  the  systems 
being  forecast. 

In  order  to  make  a  forecast,  either  *f  or  r  must  be  eliminated 

k 


from  equation  (3).   'Arnason  has  used  a  geostrophic  relation  to  eliminate 
Y  ,  obtaining  equation  (U)  which  involves  ¥     alone.  To  apply  this  equa- 
tion, height  fields  must  first  be  "balanced"  to  obtain  the  stream  function 
and  then,  at  the  conclusion  of  the  forecast,  inverted  to  obtain  the  fore- 
cast height  field.  Another  approach  would  be  to  compute  the  vorticity 
and  winds  directly  from  the  height  values  using  the  geostrophic  rela- 
tionship, obtaining  the  forecast  equations 

(Vl-£l\rttlV'V7-£lV'Vi  =  0  (6) 

\     JflVJf     %        I   W 


3.   Procedure 

As  mentioned  above,  in  order  to  apply  the  stratified  model  to  actual 
forecasting,  values  must  be  chosen  for  the  parameters  V   and  g'H.  For 
the  purposes  of  the  tests  presented  here,  it  was  decided  to  use  for  v 
a  mean  wind  computed  from  a  height  field  obtained  by  first  heavily 
smoothing  the  500-mb  and  700-mb  heights,  then  computing  an  intermediate 
field  by  the  equation 

*  ^^t^t  '25  f      •  (7) 

SOO  TOO 

The  wind  obtained  by  this  method  is  thus  a  space-mean  wind  for  some 
intermediate  level  between  500  mb  and  700  mb.  For  the  first  six  one- 
hour  time  steps,  a  V  computed  from  the  initial  data  was  used;  there- 
after a  new  V   was  computed  from  forecast  fields  every  sixth  time  step. 

In  the  choice  of  the  parameter  g'H,  consideration  was  given  to  the 
similarity  between  the 'forecast  equation  (U)  and  the  forecast  equation 
for  the  "Helmholtz"  barotropic 

Since  it  was  known  that  forecasts  made  with  this  model  at  5>00  mb  are  of 
consistently  good  quality,  it  was  decided  to  choose  g'H  such  that  the 
Helmholtz  term  for  the  stratified  model  would  be  the  same  as  the  Helm- 
holtz term  of  the  barotropic.  Thus  in  the  Helmholtz  term,  g'H  was 
chosen  to  be 

g'H  ■  .gjkl  (?) 

where  k  is  a  "tuning  constant"  with  a  value  of  four,  and  H^  is  the  stand- 
ard height  of  the  500-mb  surface. 

In  order  that  the  model  be  used  at  levels  other  than  £00  mb,  it  was 
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decided  to  include  a  second  "tuning  constant"  in  the  final  divergence 
term.  By  varying  this  second  constant  the  proper  amount  of  divergence 
may  be  introduced  at  each  level  to  produce  the  best  forecast,,  With 
these  values  introduced,  the  forecast  equation  for  the  stream  function 
case  becomes 

(  V'-H-f  7^  XI  +V-V7  -  t- k.P-7.  V'Vr°  O  do) 

v      tu*  '  J.t  '     ThT 

In  the  same  manner,  starting  with  equation  (6)  a  similar  equation  was 
obtained  for  the  geostrophic  version  of  the  model: 

/  Vl-±t£\  12  +LV-V7-  tkJ-2  V-V2  =  o  (n) 

The  model,  in  the  form  of  finite  difference  equations  obtained  from 
equations  (10)  and  (ll),  was  run  to  l|8-hour  forecasts  for  850  mb,  700  mb 
and  5>00  mb  for  several  randomly  selected  days.  Various  values  of  k-,  were 
used  at  each  level  for  both  the  geostrophic  and  stream-f unction  cases. 
In  addition,  for  the  6  May  1962  case  runs  were  made  on  1000-mb  data  in 
order  to  determine  the  feasibility  of  using  the  model  at  that  level. 


k*       Computational  Methods 

The  test  computations  of  the  stratified  model  described  in  this 
paper  were  made  on  a  Control  Data  Corporation  I60I4.  computer  using  a 
program  constructed  of  meteorological  subroutines  prepared  by  the  Fleet 
Numerical  Weather  Facility,  Monterey.  The  grid  used  in  the  numerical 
solution  of  the  prediction  equation  was  the  standard  octagonal  grid  used 
by  both  the  Joint  Numerical  Weather  Prediction  Unit  at  Suitland,  Md.  and 
the  Fleet  Numerical  Weather  Facility. 

Time  differencing  was  accomplished  in  one-hour  time  steps  using  the 
central  difference  formula 


♦m  *+*-,  + a  ( M)r  ,  <"> 


r 

where  A^^  is  the  computed  one-hour  change  in  geopotential.  At  the 


initial  time,  T=  O  ,  and  all  integral  multiples  of  twelve,  the  forward 
difference  formula 

tYt,  =*Y+  U*f)r  (1 

was  used. 

In  cases  where  a  stream  function  wind  was  used,  the  stream  values 
were  obtained  from  initial  height  data  through  the  solution  of  the  bal- 
ance equation  as  described  by  'Arnason  [I4J.  All  initial  height  fields 
and  fields  required  for  forecast  verification  were  obtained  from  files 
of  operational  analyses  at  the  Fleet  Numerical  Weather  Facility. 


5.   Results 

Verification  scores  for  each  of  the  cases  tested  were  determined 
by  a  FNWF  program  using  the  equations 

Pillow  =     £i  (A- B)n  (1^) 

N 


RMSE  -   I  ^,[(A-B)n-Pi)U^T  (15) 

V  N 

Here  A  represents  the  forecast  value,  B,  the  observed  value  and  N,  the 
total  number  of  grid  points,  1977. 

As  can  be  seen  from  the  equations,  the  pillow  is  a  measure  of  the 
error  in  the  forecast  of  the  mean  height  while  the  RMSE  (root  mean  square 
error)  is  a  measure  of  the  deviation  of  errors  around  the  mean  error. 

Figures  1,  2  and  3  show  plots  of  RMSE  as  a  function  of  k-j_  for  the 
various  cases  tested.  As  can  be  seen,  the  value  of  k-^  necessary  to  pro- 
duce the  best  forecast  (minimum  RMSE)  varies  somewhat  from  case  to  case 
at  each  level.  However,  the  curves  are  relatively  flat  in  the  vicinity 
of  the  minima  so  that  by  selecting  a  reasonable  value  of  k-, ,  the  quali- 
ty of  forecasts  from  one  case  to  another  will  vary  only  slightly  from 
the  optimum.  Table  1  shows  mean  values  of  RMSE  for  all  cases  obtained 
with  various  values  of  k,  at  each  level.  A  plot  of  these  values  is 
shown  in  figure  U.  From  this  plot  of  mean  error  versus  k-, ,  values  of 
k^  can  be  selected  which  will  give  the  minimum  total  error.  Table  2 
gives  a  tabulation  of  the  error  which  would  occur  in  each  case  tested  if 
k^  at  each  level  were  selected  as  the  value  which  produces  the  smallest 
mean  error. 

Forecasts  for  850  mb  made  from  initial  data  for  20  January  1962  are 
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Table  1.  Mean  Verification  Scores  of  Forecasts  With  Various  Values  of  k-^ 


STREAM-FUNCTION  FORECASTS 

kl 

RMSE  (ft) 

kl 

RMSE  (ft) 

kl 

RMSE  (ft) 

kl 

RMSE  (ft) 

CO 

850 

2 

112 

3 

107 

h 

108 

5 

nil 

3 
o 

700 

1 

133 

2 

12U 

3 

123 

h 

128 

C\J 

500 

0 

lli8 

i 

2 

1U7 

1 

1U7 

2 

150 

CO 

850 

2 

17U 

3 

166 

h 

166 

5 

173 

g 

700 

1 

189 

2 

17U 

3 

169 

h 

173 

CO 

1  -J- 

5oo 

0 

219 

i 

2 

215 

1 

2lh 

2 

218 

GE0STR0PHIC  FORECASTS 

kl 

RMSE  (ft) 

kl 

RMSE  (ft) 

kl 

RMSE  (ft) 

kl 

RMSE  (ft) 

CO 

850 

2 

130 

3 

123 

h 

125 

5 

13U 

3 

- 

700 

1 

137 

2 

127 

3 

126 

h 

131 

i  -J 

C\J 

5oo 

o 

163 

i 

2 

161 

1 

162 

2 

169 

CO 

850 

2 

198 

3 

192 

h 

193 

5 

200 

3 

o 

700 

1 

2 

203 

3 

195 

h 

198 

i     CO 

1  -=J 

5oo 

0 

— 

i 

2 

271 

1 

269 

2 

277 
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Table  2.  Verification  Scores  For  Each  Day  Using  Optimum  Value  of  k]_ 


850  mb 
RMSE  (ft) 

kx=3.5 

700  mb 
RMSE  (ft) 

kx-2.75 

500  mb 
RMSE   (ft) 

k^O.9 

1*  November  1961 
(2k   hr) 

111 

116 

131* 

20  January  1962 
(21*  hr) 

128 

128 

178 

6  May  1962 

(21*  hr) 

87 

88 

127 

k  November  1961 
(1*8  hr) 

173 

185 

226 

20  January  1962 
(1*8  hr) 

198 

185 

250 

6  May  1962 
(1*8  hr) 

11*3 

11*0 

168 
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shown  in  figures  8  through  13.  Initial  charts  and  charts  for  the  veri- 
fying time  of  each  forecast  are  shown  in  figures  5>  through  7.  Comparing 
forecasts  for  the  various  values  of  k^  at  a  given  level  shows  that  the 
forecast  speed  of  synoptic  waves  increases  as  the  value  of  k^  is  increased 
(thus  decreasing  the  value  of  gfH  and  increasing  the  magnitude  of  diver- 
gence). This  is  in  accordance  with  the  results  of  equation  (5).  Further, 
it  can  be  seen  that  there  is  very  little  change  in  intensity  of  the  vari- 
ous systems  as  k-j_  is  varied.  This  demonstrates  the  point  that  the  addi- 
tional divergence  term  of  the  stratified  model  results  in  movement  rather 
than  development. 

As  can  be  seen  from  equations  (10)  and  (ll)  a  choice  of  k]_=0  would 
reduce  the  forecast  to  the  barotropic.  An  examination  of  figure  k   re- 
veals that  the  optimum  value  of  k]_  at  £00  mb  is  close  to  zero,  but  that 
the  stratified  model  does  show  a  slight  improvement  over  the  barotropic 
with  values  of  k-^  between  zero  and  one.  This  is  as  expected,  since  £00 
mb  is  generally  accepted  to  be  at  or  near  the  level  of  minimum  divergence, 
and  thus  the  contribution  of  divergence  to  the  change  in  the  vorticity 
pattern  at  this  level  should  be  small. 

At  levels  other  than  5>00  mb,  where  the  barotropic  model  is  less  ade- 
quate, the  values  of  k-^  for  which  the  minimum  RMSE  is  reached  are  greater. 
For  the  700-mb  level  optimum  values  for  the  cases  tested  range  from  k-^s2 
to  kj=Us   while  at  8^0  mb  they  range  from  3  to  5.  These  results  can  be 
accounted  for  by  the  vertical  distribution  of  velocity  divergence  which 
generally  increases  with  distance  from  the  5>00-mb  surface. 

By  use  of  equation  (1)  and  the  chosen  values  for  g'H  given  by  equa- 
tion (9),  an  expression  for  the  divergence  at  a  given  time  may  be  ob- 
tained: 
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dwv  =  q-f     /!£  +  k,  v- Vf>\  (16) 

In  figures  lit,  15  and  16  divergence  charts  at  85>0,  700  and  £00  mb  as  com- 
puted from  equation  (16)  for  initial  data  of  20  January  1962  are  shown. 
These  charts  show  clearly  the  vertical  distribution  of  divergence,  with 
values  as  large  as  +  5  x  10~°  sec"  at  8^0  mb,  +  3  x  10"°  sec   at  700 
mb  and  values  very  nearly  zero  at  5>00  mb. 

An  examination  of  the  divergence  fields  relative  to  their  correspond- 
ing height  fields  reveals  that  lines  of  zero  divergence  very  closely  fol- 
low the  short-wave  trough  and  ridge  lines  of  the  corresponding  height 
fields.  Between  the  short-wave  trough  and  ridge  lines  are  located  the 
centers  of  maximum  convergence  and  divergence,  with  divergence  ahead  of 
ridges  and  convergence  ahead  of  troughs,  as  might  be  expected. 

For  the  case  of  6  May  1962,  1000-mb  forecasts  were  made.  From 
figure  3  it  can  be  seen  that  the  curves  of  RMSE  versus  k-^  for  1000  mb 
are  almost  exactly  parallel  to  the  curves  for  8^0  mb  with  an  RMSE  about 
ten  feet  greater  at  1000  mb  than  at  8^0  mb.  The  fact  that  the  curves 
are  nearly  parallel  suggests  that  the  increased  error  at  1000  mb  is  due 
to  a  greater  change  in  the  intensity  of  systems  at  this  level  than  at 
850  mb.  The  movement  of  established  systems  which  appear  on  both  initial 
charts  is  nearly  the  same  at  both  levels;  thus  k-^  which  controls  the  mag- 
nitude of  the  additional  movement  term  should  reasonably  be  chosen  the 
same  for  the  two  levels.  The  verifying  charts  for  the  forecasts  in 
question,  figures  21,  22,  25  and  26,  show  this  to  be  the  case.  Forecasts 
at  both  levels  have  the  major  systems  placed  very  close  to  their  actual 
positions  at  verifying  time,  but  the  verifying  charts  at  1000  mb  differ 
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more  in  pattern  due  to  the  variation  in  intensitjr  change  of  different 
systems  over  the  forecast  period. 

In  addition  to  the  forecasts  made  for  the  three  days  using  stream 
functions  obtained  from  the  initial  height  fields,  for  two  of  the  days 
forecasts  were  made  directly  from  the  height  fields  using  equation  (ll). 
As  can  be  seen  from  figures  27  through  30,  the  geostrophic  forecasts 
thus  obtained  showed  great  similarity  to  the  forecasts  made  from  the 
stream  function.  However,  as  expected,  the  geostrophic  version  appears 
to  suffer  from  the  major  error  common  to  geostrophic  forecasts  -  the 
"blow-up"  of  low- latitude  high-pressure  areas.  As  can  be  seen,  the  I4.8— 
hour  geostrophic  forecast  increases  the  central  height  of  the  5>00-mb 
high  off  the  east  coast  of  the  United  States  to  over  19,500  feet,  while 
the  stream-function  version  keeps  the  central  height  to  a  more  reasonable 
value.  The  treatment  of  the  high  off  the  west  coast  of  the  United  States 
is  similar. 


1U 


6.   Conclusions 

The  stratified  model  as  proposed  by  'Arnason  appears,  on  the  basis 

of  the  test  cases  presented  here,  to  produce  reasonable  forecasts  for 
levels  at  and  beloitf  £00  mb.  The  stream-function  version  of  the  strati- 
fied model  produces  500-mb  forecasts  for  2k   and  i|8  hours  which  are  in 
all  cases  very  slightly  better  than  those  produced  by  the  barotropic 
model.  At  850  mb  and  700  mb  the  model  produces  forecasts  of  system 
movement  which  are  very  realistic.  The  geostrophic  version  of  the  model 
suffers  from  the  "blow-up"  of  low- latitude  high  pressure  cells,  as  is 
common  with  most  geostrophic  models. 

Optimum  values  of  the  multiplying  factor  of  the  divergence  term 
k;j_  at  8^0,  700  and  500  mb  have  been  found  to  be  3.5,  2.75  and  0.9 
respectively.  Although  the  number  of  cases  tested  was  not  large,  the 
values  determined  here  should  be  sufficiently  accurate  to  provide  a 
basis  for  the  long-term  testing  necessary  before  any  model  becomes  suit- 
able for  operational  use. 
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Figure  £.  8£0-mb  analysis,  20  January  1962 


20 


Figure  6.  8^0-mb  analysis,  21  January  1962 
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Figure  7.  850-mb  analysis,  22  January  1962 
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Figure  8.  2U-hour  forecast  for  8^0  mb  from  20  January  1962,  k-^0 
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Figure  9.  2U-hour  forecast  for  850  mb  from  20  January  1962,  k  =3 
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Figure  10.  2l;-hour  forecast  for  8£0  mb  from  20  January  1962,  k  =$ 
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Figure  11.  U8-hour  forecast  for  350  mb  from  20  January  1962,  k,=0 
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Figure  12.  I|8-hour  forecast  for  8£0  mb  fr9m  20  January  1962,  k  =3 
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Figure  13.  US-hour  forecast  for  8£0  mb  from  20  January  1962,  k  =£ 
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Figure  llw  850-mb  divergence,  20  January  1962 
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-6    -1 
contour  interval  =  2  x  10   sec 


Figure  15.  700-mb  divergence,  20  January  1962 
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-6   -1 
contour  interval  =  2  x  10   sec 


Figure  16.  500-mb  divergence,  20  January  1962 
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Figure  17.  1000-mb  analysis,  6  May  1962 
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Figure  18.  8£0-mb  analysis,  6  24ay  1962 


33 


'■■■'.    '    '.'•■?' 


Figure  19.  2l*-hour  forecast  for  1000  mb  from  6  May  1962 
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Figure  20.  2U-hour  forecast  for  8^0  mb  from  6  May  1962 
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Figure  ZL.     1000-mb  analysis,  7  May  1962 
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Figure  22.  850-rab  analysis,  7  May  1962 
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Figure  23.  Ii8-hcmr  forecast  for  1000  mb  from  6  May  1962 
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Figure  2l*.  l|8-ho-ur  forecast  for  8^0  mb  from  6  May  1962 
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Figure  25.  1000-mb  analysis,  8  May  1962 
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Figure  26.  850-mb  analysis,  8  May  1962 
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Figure  27.  500-mb  analysis,  h   November  l°6l 
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Figure  28.  U8-hour  forecast  for  £00  mD  from  h   November  1961 
(stream  function) 
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Figure  29.     1+3-hour  forecast  for  500  mb  from  k  November  19  61 
(geostrophic) 
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Figure  30.  500-mb  analysis,  6  November  1961 
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